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We revisit the extraction of aa(A/^) from the QCD perturbative corrections to the hadronic r 
branching ratio, using an improved fixed-order perturbation theory based on the exphcit summation 
of all renormalization-group accessible logarithms, proposed some time ago in the literature. In this 
approach, the powers of the coupling in the expansion of the QCD Adler function are multiplied 
by a set of functions D„, which depend themselves on the coupling and can be written in a closed 
form by iteratively solving a sequence of differential equations. We find that the new expansion has 
an improved behaviour in the complex energy plane compared to that of the standard fixed-order 
perturbation theory (FOPT), and is similar but not identical to the contour-improved perturbation 
theory (CIPT). With five terms in the perturbative expansion we obtain in the MS scheme aa{M'^) = 
0.338 ±0.010. 
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I. INTRODUCTION 

The non-strange hadronic decays of the r lepton pro- 
vide one of the most precise determination of the strong 
coupling Us- The recent calculation of the Adler func- 
tion to four loops the same order to which the /3- 
function of the renormalization-group (RG) equation is 
known [2, Jsj, renewed interest in the determination of 
ctcfMj) |4]-|20l|. The intriguing remark Q that the inclu- 
sion of a higher order term increased, instead of reducing, 
the theoretical error on the resulting as(M^) stimulated 
many investigations aimed at understanding this fact. 

The basic procedure involves the analytic continuation 
of the Adler function (the logarithmic derivative of the 
massless QGD polarization function) in the complex en- 
ergy plane, where it can be calculated by the Opera- 
tor Product Expansion (OPE). The contribution of the 
higher dimensional terms ("power corrections") in the 
OPE was evaluated and found to be quite small [l-Q. 
Recently, also deviations of the true polarization function 
from the OPE description, especially near the timclike 
axis (the so-called "duality violation terms" ) were evalu- 
ated in a model and were shown to have small effects on 
the coupling determination j20| . Thus, the problem has 
to do mainly with the pure perturbative part of QCD, 
more exactly with the notorious ambiguities of the trun- 
cated perturbative series. 

There are two competing versions of perturbation 
theory, the standard fixed-order perturbation theory 
(FOPT) and the RG improved, which in this context 
is also known as contour-improved perturbation theory 
(CIPT). Their predictions differ by about 0.02, which 
is at present the main part of the theoretical error on 
as{M^) i IE [11 [11. 

The investigation of the perturbative series of QCD 
in the context of the uncertainty in the extraction of 
is of such great importance that its theoretical aspects 
have been studied by several authors and various alter- 
native approaches have been proposed. They include in 
general additional information about the scries beyond 



the truncation order, known either from specific classes 
of Feynman diagrams or from RG invariance. Thus, a 
reordering of the standard contour-improved approach 
exploiting RG invariance was proposed in and a de- 
tailed analysis of the errors of various expansions has 
been performed in Q. 

A more radical modification was investigated in [3, [isl . 
[l3 | , where the available knowledge on the large-order be- 
havior of the perturbative coefficients was exploited with 
mathematical techniques of accelerating the series con- 
vergence by means of conformal mappings j2l| - [2^ . This 
led to a modified expansion in terms of a new set of func- 
tions, which have the advantage of sharing the known 
singularities of the expanded correlator in the coupling 
and the Borel complex planes. As argued in this ex- 
pansion is particularly suitable in the contour-improved 
version, since it make a summation of both the running 
coupling and of the Feynman coefficients of the Adler 
function. Detailed numerical studies 0, [l!] proved the 
good convergence properties of the latter expansion for a 
large class of physical models which simulate the known 
properties of of the Adler function. 

In the light of the above, any fresh attempt to improve 
the understanding of the properties of the perturbative 
expansion in the complex energy plane and the origin of 
the discrepancy in the coupling predictions would be wel- 
come. With this motivation, wc consider in the present 
paper a RG improved expansion proposed in l25l . [26j . 
using a procedure originally advocated in [27l429j . The 
method is a generalization of the leading logarithms sum- 
mation, in which terms in powers of the coupling constant 
and logarithms are regrouped, so that for a given order, 
the new expansion includes every term in the perturba- 
tive scries that can be calculated using the RG invariance. 
The method was applied to several correlators and ob- 
servables, for instance the inclusive decays of the 6-quark 
[25| and the hadronic cross section in e+e~ collisions [26| . 
where its main merit was proved to be a substantial re- 
duction in sensitivity to the renormalization scale. In 
the present work we investigate the new expansion for 
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the QCD Adlcr function in the complex energy plane 
and the determination of as from r hadronic decays. To 
our knowledge, this problem was not investigated in full 
generality up to now.^ 

The scheme of this paper is as follows: for completeness 
we briefly review in Sec. |TT] the perturbative expansion 
of the Adler function and its connection to the hadronic 
decay width of r. In Sec. IIIIl following Rcf. [2^, we re- 
view the derivation of the new RG improved expansion of 
the Adler function and give the corresponding expansion 
functions calculated to four loops. For further applica- 
tions of the method it is useful to know also the higher 
expansion functions, which we have calculated in an an- 
alytic closed form by iteratively solving the relevant dif- 
ferential equations. As the general expressions are rather 
lengthy, we give in the Appendix simpler forms of the ex- 
pansion functions up to n = 10 obtained by inserting the 
numerical values of the known perturbative coefficients 
of both the Adlcr and /3-functions to four loops. The ex- 
pressions are written in terms of the coefficients beyond 
four loops, which are not yet available from explicit cal- 
culations and are left arbitrary. In Sec. llVl we investigate 
the properties of the new expansion in the complex en- 
ergy plane and compare it with the standard FOPT and 
CIPT, using in particular a physical model for the Adlcr 
function proposed in Q. In Sec. |V] we report a new 
determination of as{M^) based on the FO expansion im- 
proved by RG summation discussed in this paper. Sec. 
I VII summarizes our results and presents some conclusions. 



II. ADLER FUNCTION IN PERTURBATIVE 
QCD 

The Adlcr function plays a crucial role in the determi- 
nation of as{M^) from hadronic r decays. The method 
is discussed in the seminal paper [30l | and is reviewed in 
several recent articles [1, [^, [l^, [l^ . For completeness we 
give below a few details. 

The inclusive character of the total r hadronic width 
makes possible an accurate calculation of the ratio 

_ r[T^ — > i^T-hadrons] 

lir — 



Of interest is the Cabbibo allowed component which pro- 
ceeds either through a vector or an axial vector current, 
since in this case the power corrections are particularly 
suppressed. On the theoretical side, i?r can be expressed 
in the form 



D>2 



, (1) 



^ The RG summation discussed in [271.12811 has been apphed to the 
extraction to as from t decays in l29l| , but only using the pertur- 
bation series to NNLO treated with Borel summation methods. 



where Nc = 3 is the number of quark colours, S'ew 
and are clectrowcak corrections, (5^°^ is the domi- 
nant perturbative QCD correction, and the (5^^^' denote 
quark mass and higher D-dimensional operator correc- 
tions (condensate contributions) arising in the OPE. 

Unitarity implies that the inclusive hadronic decay rate 
can be written as a weighted integral along the timelike 
axis of the spectral function of the polarization function 
n(^+")(s), where the superscript denotes the spin. As 
shown in [soj . the analytic properties of the polariza- 
tion function and the Cauchy theorem allow one to write 
equivalently this quantity as an integral along a contour 
in the complex s-plane (chosen for convenience to be the 
circle \s\ = M^). After an integration by parts, in our 
notation the quantity of interest 5^°^ is expressed as the 
following contour integral: 
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where the reduced function D{s) = D^^^'^^{s) — 1 is ob- 
tained by subtracting the dominant term from the Adler 
function, i.e. the logarithmic derivative of the polariza- 
tion function, D(i+o)(s) ee -s dn(i+o)(s)/ds 

The function D{s) depends only on the energy variable 
s, but for its pure perturbative part Dpcrt appearing in 
([2]) we emphasized also the formal dependence on the 
renormalization scale entering through the coupling 
and the standard perturbative logarithms. Specifically, 
we shall denote by 
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In the so-called "fixed-order perturbation theory", one 
chooses a fixed scale = and the expansion of D 
reads 



£>FOPT(a, L) = 



n=l fc=l 



fc-1 



(4) 



In the expansion above, the leading coefficients c„_i are 
calculated from Feynman diagrams. The known coeffi- 
cients Cn,i are (see [l| and references therein): 

ci.i = 1, C2,i = 1.640, C3,i = 6.371, C4.1 = 49.076, (5) 

and several estimates for the next coefficient 05^1 were 
made recently [6, 18, 19] . The remaining coefficients Cn,k 
for fc > 1 are determined from RG invariance and in- 
volve the coefficients Pj appearing in the perturbation 
expansion of the RG /3-function 
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The /3-function was calculated to four loops in the MS- 
renormalization scheme, the known coefficients being (see 
for the calculation of and earlier references): 

/3o = 9/4, /3i = 4, /32 = 10.0599, /Jg = 47.228. (7) 
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As remarked in [3l|, due to the large imaginary parts 
of the logarithms of —s/M"^ along the circle \s\ = M^, 
the series (|4]) is badly behaved especially near the time- 
like axis. This mandates one to search for expan- 
sions that would be better behaved and would exhibit a 
smaller renormalization-scale dependence. The "contour- 
improved perturbation theory" [3l|, |32| is based on the 
RG improved expansion, defined by the choice = —s, 
when (U) reduces to 
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The main improvement comes from the treatment of the 
running coupling as(— s), which is determined by solving 
the RG equation ^ numerically in an iterative way along 
the circle, starting with the input value as(M^) at s = 
-Ml 

The expansions (0]) and ^ coincide formally as long as 
all the terms in the series are retained (we ignore in this 
discussion the fact that the coefficients c„_i are known 
to increase as n! and the series are actually divergent). 
However, since the expansion coefficients are known only 
up to a finite and not so large order, the series have to 
be truncated at some order n < N. Then the expansions 
differ by terms of order a^~^^, which may be substantial 
due to the relatively large value of the coupling at the 
low scale set by the mass of the r. Therefore, the expan- 
sions lead to different values for S^'^\ this being the main 
source of error in the determination of as(A/^) from the 
hadronic T-dccays. 



III. RENORMALIZATION GROUP 
SUMMATION 

As suggested in [H, [2^ , the FO expansion dU of the 
reduced Adler function can be written equivalently as 



(9) 



where the functions Dn{u), depending on a single vari- 
able u = aL, are defined as 



Dn{u) 
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As seen from the definition, the first function Di sums all 
the leading logarithms, the second function D2 sums the 
next-to-leading logarithms, and so on. Thus, the sugges- 
tion was to effectively make a summation by collecting 
the aggregate coefficients of the leading logarithms mul- 
tiplied by fixed powers of the coupling constant. The 
attractive feature pointed out in [251 |26[, is that these 
functions can be obtained in a closed analytical form. 
We sketch below the derivation, which is based on RG 
invariance. 

The Adler function defined by ^ , calculated in a fixed 
renormalization scheme, is scale independent and satisfies 



the RG equation 

M'^{5FOPT(a,i)} -0, (If) 
which can be written equivalently as 
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/3(a) — = 0. 



da 



dL 



(12) 



Using in this relation the expansion ^ yields the follow- 
ing equation: 
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n=l k=2 

(/3oa^ + /3ia^ + ^2a'* + • • • + /3;a'+^ + . . .) 

00 n 

x^^7ifcc„,fca"-iL'^-i. (13) 
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By extracting the aggregate coefficient of a^L" p one 
obtains the recursion formula {n > p) 

p-2 

0= {n-p + 2)c„,„_p+2 + ^(?^ - i - l)PtCn-t-l,n-p+l- 

1=0 

(14) 

These relations are well known, and in particular for n < 
4 they coincide with the relations given in eq. (2.11) of 

i- 

Multiplying both sides of (O by {n-p + 1)u"~p and 
summing from ri = p to 00, we obtain a set of first-order 
linear differential equation for the functions defined in 
([To)), written as 



p-2 



p-2 



du ^ du 



Setting now n = p — 1 we write this set as 

n-l 

R„ ( i„ 

du 



(15) 



^ + 'f2l3e(u^ + n-£]D,,.e = 0, (16) 

du -"^ — ' \ nil I 

(=0 

for ri > 1, with the initial conditions I?„(0) ~ c„^i which 
follow from ((TUj). 

The solution of the system ([T6[) can be found itera- 
tively in an analytical closed form. It turns out that the 
solutions Dn(u) depend on u only through the variable 
w = 1 -\- (3ou. The expressions of Dn{u) for 71 = 1, 2, 3, 4, 
written in terms of this variable, and the coefficients c„_i 
and j3k are: 

Cl,l 



Di{u) = — '-, w = 1 -I- /3ou, 

w 
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'0 /^o 

(-/33 + /32/33)ci,i 2(-^2+/3o /32)C2, 



C4,l 



/3i(-2/32ci,i + 3/3o/32Ci,i+2/3o/3iC2,i+3/3^C34)lnw _4 ^ ^^(5^^^^^^ _^ g^^^^^^) ^^^2 ^ ^g^d^iln^w _4 
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In [251 126| similar differential equations were solved for 
n < 4 for several observables, including the cross section 
of e+e^ annihilation into hadrons, whose expansion in 
QCD is related to the expansion of the Adler function 
in which we are interested. The functions D„{u) given 
above coincide actually with those calculated in [2g] . For 
the applications made in this work and possible further 
studies, we have derived the expressions of -D„ up to n = 
10. The solutions depend on the coefficients c„^i and 
the coefficients (3k of the expansion ^ of the /3-function. 
For consistency, to each Feynman diagram order n we use 
the expansion of the /3-function to the same order. The 
complete expressions are rather lengthy. They simplify 
considerably if we insert the known numerical values of 
the coefficients c„,i for n < 4 given in (O, and of the 
coefficients f3k for fc < 3 given in ([7]). The corresponding 
expressions, which depend on the arbitrary coefficients 
c„^i for 5 < n < 10, and /3k for 4 < fc < 9, are listed in 
Appendix. 

We shall use in what follows the truncated FOPT im- 
proved by renormalization-group summation (RGS) writ- 
ten as 
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IV. DISCUSSION 

In this section we shall investigate the properties of the 
new expansion (jl9p in the complex s-plane, along the cir- 
cle s = exTp{i6). For comparison, we plot in Figs. [T][3] 
the modulus of each successive term of order n < 5 of 
the standard FOPT expansion the CIPT expansion 
^ and the RGS improved FOPT expansion (jl9p , respec- 
tively. For convenience, we have taken as(M^) = 0.34. 
For n — 5 we have used the expression of given in the 
appendix, with the estimate cs^i ~ 283 from ITgj and 
setting (34 = 0. 

From Fig. [1] it is seen that the higher order terms 
are large near the timelike axis {6 = 0). This shows the 
slow convergence of the standard FOPT in this region. 
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FIG. 1: Modulus of the perturbative terms of the standard 
FO expansion Q along the circle s = exp{i9). 



where the logarithms defined in ([3|) acquire large imagi- 
nary parts. As discussed in [3lj, the reason is the poor 
convergence, especially near the timelike axis, of the ex- 
pansion of as{—s) in powers of as{M^), which is used in 
passing from the renormalization-group improved expan- 
sion to the fixed-order expansion ((T5)) . In contrast. 
Fig. [5] shows that in CIPT the higher terms are much 
smaller, i.e. the expansion has a good convergence along 
the whole circle. As seen from Fig. [H the RGS improved 
FOPT expansion ([T^ has a behavior similar to that of 
CIPT: the series is stable along the circle and the higher 
order terms are very small. Thus, although it depends 
explicitly only on the coupling at a fixed scale, the expan- 
sion shares the good qualities of the CI expansion 
along the circle, as seen from Fig. SJ where we simulta- 
neously plot the first three terms for the two expansions. 

In order to see the difference between CIPT and the 
FOPT improved by RGS, it is useful to look at the lead- 
ing term, with n ~ I. In the CI expansion (|5|) this term 
is ci_iQ!s(— s)/7r, where the coupling is calculated as the 
numerical solution of the RG equation ([6]) , keeping four 
terms in the expansion of the /?-function. On the other 
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FIG. 2: Modulus of the perturbative terms of the CI expan- 
sion ([8} along the circle s = exp{iO). 



FIG. 4: Comparison of the CI expansion (O and the improved 
FO expansion p9p along the circle s = exp(iS). 
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FIG. 3: Modulus of the perturbative terms of the improved 
FO expansion ()19|l along the circle s = exp{i6). 



FIG. 5: Adler function expansions ((4]), ^ and ()19|) . summed 
up to the order N = 5, along the circle s — exp{i6). 



hand, using (fT9|) and (fT7|) wc write the leading term of the 
RGS improved expansion as ci,ia/(l + /3oa hi(— s/Af.^)) 
where a = as(M^)/7r. This is actually the exact solution 
of the RG equation (O to one loop, written in terms of 
the input as{M'^). The similar behavior of the curves 
corresponding to n = 1 in Fig. 2] shows that the effect of 
the higher order terms in the expansion of the /? function 
is small. Moreover, the smallness of the next terms of the 
expansion (fT9|) proves that the summation of the leading 
logarithms is very efficient also to higher orders. 

Fig. Oshows the behaviour along the circle of the Adler 
function given by the first N ^ 5 terms in the expansions 
(|4|), ([8]) and (fT9|) . respectively. The new FO expansion 
improved by RGS is very similar to the CI expansion, as 
expected from the previous figures. 

By inserting the FOPT, CIPT and RGS improved 
FOPT expansions (|4]), ^ ([T9)) . respectively, truncated 
at some TV, into the definition ^ of S^^\ we obtain 

the corresponding values denoted as <^poPT' '^cipt ^^'^ 



(5jpQprp respectively. In Table |I] we list these values for 
various truncation orders < 5, using in the calcula- 
tion the standard value as{M^) = 0.34. As remarked 
already, CIPT shows a faster convergence compared to 
the standard FOPT. To order A^ = 4, the difference be- 
tween FOPT and CIPT is 0.0215, which, as remarked, is 
the dominant theoretical uncertainty in the extraction of 
as from the hadronic r decay rate. On the other hand. 
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= 1 


0.1082 


0.1479 


0.1455 
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= 2 


0.1691 


0.1776 


0.1797 
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= 3 


0.2025 


0.1898 


0.1931 


N 


= 4 


0.2199 


0.1984 


0.2024 


N 


= 5 


0.2287 


0.2022 


0.2056 



TABLE I: Predictions of by the standard FOPT, CIPT 
and the RGS improved FOPT, for various truncation orders 
N. 
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for TV = 4. the difFerence between the results of the RGS 
improved FOPT and the standard FOPT is 0.01754, and 
the difFerence from the RGS improved FOPT and CIPT 
is 0.0039, which confirms that the new expansion gives 
resuhs close to those of the CIPT. For iV = 5, using the 
estimate C5,i = 283 from Q, we find that the RGS im- 
proved FOPT differs from FOPT by 0.0232, and from 
CIPT by 0.0035. 

It is of interest to see whether this behaviour is pre- 
served to higher orders. To this end we consider a class 
of physical models of the Adler function used for testing 
various expansions in 0, 0, [H, [13, [I!] . 

In particular, we consider the model proposed in Q, 
where the Adler function is defined in terms of its Borel 
transform B{u) by the principal value prescription 



5(s) = -ipv 

Po 



(20) 



where the function B{u) is expressed in terms of a few 
ultraviolet (UV) and infrared (IR) renormalons 

Bb.^{u) = SU^H+B^'^H+Bf (21) 

In ml these terms were written as 



(p — u)^p 
(p + u)"'p 



l + bi{p-u) + 



l + bi{p + u) + 



(22) 



where most of the parameters were obtained by impos- 
ing RG invariance at four loops. Finally, the free param- 
eters of the model, i.e. the residues dY^,d^ and d]^ 
of the first renormalons and the coeficients d^'^^d^'^ of 
the polynomial in (|2ip . were fixed by the requirement of 
reproducing the perturbative coefficients Cn,i for n < 4 
from ([5]) and the estimate C5.1 = 283, and read: 



expressions of D„ given in the appendix, setting /3fe = 
for fc > 4 as in the previous similar analyses of higher 
order expansions @, 0, [H, [11] . 

The figure shows that the FOPT improved by RGS 
gives results close to the CIPT predictions at all orders 
up to iV = 10. In fact, as remarked in for this partic- 
ular model the standard FO expansion describes better 
than the CIPT the "true" function. Indeed, as seen in 
Fig. [HI up to = 10 the predictions of the CI expansion 
stay below the true result, and in fact never approach it 
(for higher truncation orders A^ all the three expansions 
start to show big oscillations, due to the divergent charac- 
ter of the series). For other physical models, discussed in 
[l3 . [TgI , the CI expansion may give better results at low 
orders. In the present discussion we are not concerned 
with the validity of such models: we have used them 
only in order to compare the RGS improved FO expan- 
sion with the standard FO and CI ones at higher orders. 
The conclusion, which is model- independent, is that the 
RGS improved fixed-order expansion gives results close 
to those obtained by the contour-improved expansion. 
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FIG. 6: Dependence of J*"' in FOPT, CIPT and RGS im- 
proved FOPT on the truncation order A'' in the Beneke & 
Jamin model [^. The gray band is the exact value obtained 
with the expressions (|20|) - (|23p . 



d^° = 0.781, = 7.66 X 10"^ (23) 

Then all the higher order coefficients c„_i are fixed and 
exhibit a factorial increase, showing that the perturbative 
series of the Adler function is divergent. We list below 
the values which we use in our analysis 

C5,i = 283, C6,i = 3275, 07,1 = 1.88 x 10^, 

cs.i = 3.88 X 10^ C9,i = 9.19 x 10^ cioa = 8.37 x 10^. 

In Fig. [6j we show the exact value of iJ'"' obtained with 
the above model, and the dependence of the truncation 
order A^ for the three expansions considered: standard 
FOPT, standard CIPT and RGS improved FOPT. As 
in the previous figures we have used as input as (7\jf^) = 
0.34. For the RGS improved FOPT we have used the 



V. DETERMINATION OF as{M?) 

In this section we shall use the RGS improved FO ex- 
pansion (jlOp for a new determination of as{M^) in the 
MS scheme. We use as input the recent phcnomenological 
value of the pure perturbative correction to the hadronic 
T width [H 

Ccn = 0-2037 ± 0.0040e.p ± 0.0037pc, (24) 

where the first error is experimental and the second re- 
flects the uncertainty of the higher order terms ( "power 
corrections") in the OPE. 

For the theoretical evaluation of (5'°-' from ([2]) we ap- 
ply the improved FO expansion truncated at A^ = 5, 
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choosing the scale as = ^M^ with ^ = 1 ± 0.5. We 
have used the functions D„ for n < 5 given in Ap- 
pendix, taking as input in the conservative estimates 
C5,i = 283 ± 283 Hi 0, [H and /34 = ± /3|//32, as in 
[j, [3 [l^ • With this input we obtained from the above 
phenomenological value of S^^'^ the prediction 

«,(Af2) = 0.3378 ± 0.0046exp ± 0.0042pc i;];^^?^(c5,i) 

+0.0005/ 1 \ ,+0.000085 / o ^ /or;-, 

In this result the first two errors are due to the corre- 
sponding uncertainties of 5^^^^^ given in (j24p . the third 
one reflects the uncertainty of the coefficient cs^i with 
the very conservative range adopted above, the fourth is 
due to scale variation, and the last one shows the effect 
of the truncation of the /3-function expansion. One may 
note the very small sensitivity of as(M^) on the variation 
of the scale, and a relatively large contribution of the un- 
certainty of the five loop coefficient C51 , a feature noticed 
also in the standard CIPT analyses [isl . [Tsj and in the CI 
expansions improved by the conformal mappings of the 
Borel plane 0,[ll. 

Combining in quadrature the errors given in (j25p . we 
write our prediction as 

a^(Af^) = 0.338 ±0.010. (26) 

We mention that for the same input the standard FOPT 
and CIPT give, respectively, 

a,(M2) = 0.320t°°J?, FOPT, 
a,(M2) = 0.342 ±0.012, CIPT. (27) 

For comparison we mention also the value as{M^) = 
0.320to;oi4j obtained recently in [l3| by using an im- 
proved CI expansion based on the analytic continuation 
in the Borel plane. It is of interest to quote also the very 
recent lattice determination [ssj . which predicts in the 
MS scheme: a^iM^) 0.339 ± 0.013, very close to our 
result (Uni)- 

VI. SUMMARY AND CONCLUSION 

In this paper we have applied the method of explicit 
summation of all RG-accessiblc logarithms proposed in 
[25I l26j to the perturbative expansion of the Adler func- 
tion relevant for the extraction of as from r hadronic 
decays. The work is motivated by the well-known dis- 
crepancy between the predictions of as{M^) from the 
standard fixed-order and RG-improved expansions. As 
this discrepancy has to do with the behaviour of the per- 
turbative expansion of the Adler function along the con- 
tour involved in the integral ^ , especially near the time- 
like axis, it was of interest to see whether a more general 
fixed-order expansion can be found, with good conver- 
gence properties along the contour. While the method 
proposed in [2^ [2^ was applied to several other observ- 
ables, its properties in the complex energy plane were not 
investigated up to now. 



The main result of the paper is that the summation of 
leading logarithms provides a systematic expansion with 
good convergence properties in the complex plane, in- 
cluding the critical region near the timelike region. By 
summing up pieces of the standard fixed-order series (jl9p 
into the functions _D„ defined in (jlOp . the new expan- 
sions (jl9p is no longer plagued by large imaginary parts 
of the logarithms, responsible for the poor convergence 
of FOPT along the contour. 

On the other hand, the results of the new expansion are 
similar to those obtained by the CI expansion ([8]) , which 
was to be expected since both implement RG invariance. 
However, the two expansions are not identical: CIPT 
uses the exact solution of the RG equation to four loops, 
found numerically by an iterative integration along the 
circle, while the new expansion involves only analytical 
closed form expressions valid along the whole integration 
contour, thereby avoiding numerical integration. As dis- 
cussed in Sec. llVi the convergence of the new expansion 
along the circle |s| = is comparable to that of CIPT. 

As mentioned earlier, several modifications of the stan- 
dard FO and CI perturbative expansion were recently 
proposed and applied to the determination of the strong 
coupling from r decays 0, [HI, [l3l • The present approach 
exploits RG invariance in a complete way, summing in 
analytical closed expressions all the terms that can be cal- 
culated to a definite Feynman diagram order. Of course, 
the truncated expansions of the different summations dif- 
fer among each other by terms of order a^~^^, which may 
be quite important at the relatively low scale relevant in 
T decays. Moreover, the actual differences depend on 
the detailed form of including known information on the 
higher order terms. Therefore, our study contributes to 
the assessment of the ambiguities of the perturbation ex- 
pansion of the Adler function in the complex plane and 
the theoretical error of as{M^). 

Using an input adopted in recent works [H, [l^, the 
new expansion (fT9|) to five loops gives the value (|26)) for 
as{M^) in MS scheme. The result is situated between 
the predictions of FOPT and CIPT given in ([?7)) . closer 
to the latter, and remarkably, is practically identical to 
the recent lattice prediction obtained in |33| . 
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Appendix: Expressions of the functions Dn 

We give the expressions of the functions £)„(«), n = 
1,2,. ..,10 in a readily readable form using the known 
numerical values of the coefficients Cn,i for 1 < n < 4 
from (P, and of Pj for < j < 3 from The higher 
coefficients c„^i for n > 5 arc arbitrary. For generality. 
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at each order n we include the higher loop coefficients Pj 
for J > 4 up to the corresponding order. 

As remarked in Sec. Illli the functions Dn{u) depend 
only on the variable 

w = l + 9/4it. (A.l) 
The explicit expressions are: 

Di{u) = w'^. (A.2) 



D2{u) =(1.64 - 1.778 \nw)w~'^. (A.3) 



I 

D3,{u) = - l.Sllw"^ + (7.682 - 8.992 lnw + 3.160 In^w) vj-^ . (A.4) 



DAu) ^ - 5.356w^^ + (-6.629 + 4.660 In w) ly-^ 

(A 5) 

+ (61.061 - 56.954 hi w + 29.596 In^ w ~ 5.619 In^ w) w~^. 

D^{u) = (20.740 - 0.148/34) w''^ + (-25.371 + 19.043 In w) w'^ 
+ (-41.986 + 43.637 In w - 12.426 In^w) w^^ 

. (A.6) 

+ (46.618 + 0.148/34 + C5,i - 535.458 In w) 

+ (255.117 In^ w - 80.143 In^ w + 9.989 In'^ w) w'^ 

Deiu) = (-8.802 + 0.395/34 - O.Ul/^s) 

+ (118.935 - 0.749/34 + (-73.7407 + 0.527/34) Inu;) 
+ (-155.498 + 169.168 In u; - 50.782 In^ w) w'^ 

(A 7) 

+ (-394.738 + 376.142 In u; - 177.243 In^ u; + 29.455 In^ w) w"^ 

+ (440.104 + 0.354/^4 + 0.111^5 + ce.i + (-1366.3 - 1.317/34 - 8.889c5,i) \nw) w'^ 

+ (2833.36 In^ w - 898.378 In^ w + 195.853 In'' w - 17.758 In^ w) u;"^ 

Dr{u) = (1.850 - 0.048^4 + 0.316^5 - 0.089/36) w""^ 

+ (-70.196 + 2.386/34 - 0.562/35 + (31.297 - 1.405/34 + 0.395/3,5) Inw) w''^ 

+ (793.632 - 4.746/34 + (-765.413 + 4.933/34) Inu; + (196.642 - 1.405/34) In^ w) w'^ 

+ (-1474.52 + 1406.51 In w - 691.764 In^ w + 120.371 In^ w) w'^ 

+ (-1007.23 - 0.971^4 - 6.553C54) w'^ 

+ (4177.48 In w - 1986.84 In^ w + 577.528 In^ w - 65.455 In'' w) w'^ 
+ (1756.47 + 3.378/34 + 0.246/35 + 0.089/36 + 6.553c5,i + C74) uj-'^ 
+ ((-7123.42 - 6.119/34 - 1.185/35 - 15.803c5,i - 10.667c6,i) In w) w"^ 
+ ((12324 + 7.023/34 + 47.407c5,i) In^ w - 11671.3 In^ w) 

+ (2743.86 In'' w - 449.389 In^ w + 31.569 In^ w) w'"^ (A.8) 
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Ds{u) = (-194.169 + 1.242/34 - 0.040/35 + 0.263^6 - 0.074/37) w'^ 
+ (-189.727 + 0.309^4 + 1.890/35 - 0.449/36) w'^ 
+ (-550.03 + 16.293/34 - 3.559/35) u;""* 

+ ((430.016 - 15.225/34 + 3.699/35) Inw + (-83.458 + 3.746/34 - 1.053/35) In^ u;"'* 
+ (7347.75 - 44.622/34 + (-7004.34 + 42.519^4) Inu;) w"^ 
+ ((3071.05 - 20.036/34) In^ u; + (-466.114 + 3.329^4) In^ w) uj-^ 
+ (-2576.43 - 3.967/34 - 26.779c5a + 15607.3 Inu;) w''^ 
+ (-7480.95 In^ w + 2263.66 In^ w - 267.492 In^ w) w'^ 
+ (-5251.36 - 4.511^4 - 0.874/35 - II.65OC54 " 7.863c64) w"'^ 
+ ((18170.5 + 10.355/34 + 69.897c5,i) Inw ~ 25812. In^ w + 8091.04 In^ w) w''^ 
+ (-1656.441x1'' u; + 139.637 In^w) w''^ 
+ (1413.96 + 35.257/34 + 2.584/35 + 0.186^6) 
+ (0.074/37 + 38.429C54 + 7.863c6.i + cg^i) w"^ 

+ ((-34522.1 - 52.922^4 - 5.167^5 - 1.106/36 - 109.64c5,i - 18.963c6,i - 12.444c7,i) \nw) w'^ 
+ ((66232.9 + 50.564/34 + 7.374^5 + 182.606c5,i + 66.370c6,i) In^ w) w'^ 
+ ((-71870.8 - 29.134^4 - 196.653c5,i) In^ +41188.6 In'^ w) w"* 

+ (-7628.07 ln^ii; + 988.189 In^w- 56.123 In^M;)^-* (A.9) 



Dg{u) = (395.544 - 10.428/34 + 0.028/3| + 1.064^5 - 0.034/36 + 0.226/37 - 0.063/38) w'^ 
+ (-462.494 - 3.971/34 + 0.022^| + O.I6O/35 + 1.565^6 - 0.375^7) w'^ 
+ ((690.377 - 4.415/34 + 0.142^5 - 0.936/36 + 0.263/37) In w) w'^ 
+ (-1810.96 + 2.864/34 + 12.852^5 - 2.848^6 - 2.220 x lO^^^^g) w^'^ 
+ ((1000.18 - 1.345/34 - 12.076/35 + 2.960/36) Inw) w'^ 
+ ((17.544 - 0.455/34 + 2.997/35 - 0.843/36) In^ w) u;""* 

+ (-5114.46 + 150.336/34 - 33.466/35 + 1.421 x l0-^^l3e - 3.553 x IQ-^^/St - 8.882 x IQ-^^Ps) w'^ 
+ ((4675.79 - 142.926/34 + 31.890/35) In u; + (-1677.31 + 60.792/34 - 15.027/35) In^ w) uj-^ 
+ ((197.827 - 8.879/34 + 2.497/35) In^ w) uj-^ 

+ (28164.8 - 98.496^4 - 0.110/3^ + 3.411 x IQ-'^^Pr, - 2.842 x lO'^^Pe) w"^ 

+ (1.421 X 10-"/37 + 3.553 x IQ-'^^Ps + 103.698c5,i - 0.741/34C5,i) w'^ 

+ ((-77765.5 + 472.226/34) In w + (36590. - 224.594/34) In^ w) w'^ 

+ ((-9928.07 + 65.284/34) In^ w + (1035.81 - 7.399/34) In'' w) w^^ 

+ (-14763 - 14.620/34 - 2.665 x 10~i^/3^ - 3.571/35) u;"^ 

+ (1.421 X 10^13/36 + 1.776 x IQ-^^/Ss - 21.844c5,i - 32.135c6,i) 

+ ((55228.2 + 42.318/^4 + 285.647c54) \nw - 96538.4 In^ w + 30623.2 In^ w) w"^ 

+ (-6511.98 In"* u; + 570.649 In^ «;) w~'' 

+ (-25449.6 - 39.014/34 - 8.882 x 10~^^^| - 3.809/35 - 0.815/36 + 7.105 x W'^^Pt) w"* 

+ (1.776 X 10"^^/38 - 80.826c5,i - 13.979c6,i - 9.174c7,i) 

+ ((97653.3 + 74.551/34 + 10.873/35 + 269.233c54 + 97.856c64) \nw) w'^ 

+ ((-158949. - 64.432/34 - 434.915c54) In^ w + 121456 In^ w) w'^ 

+ (-28116.9 In"* w + 4370.93 In^ w - 289.617 In^ w) w'^ 

+ (19040.2 + 13.328/34 + 0.060/3^ + 26.770/35 + 2.132/36 + 0.149/37) w'^ 
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+ (0.063/38 - 1.028c5,i + O.74I/34C54 + 46.115c6.i + 9.174c7.i + cg^) w'^ 
+ ((-81482.4 - 595.514/34 - 45.930/35 - 4.615/36 - 1.054^7) w'^ 
+ (-741.46c5,i - 145.547c6,i - 22.124c7a - 14.222c84) \nw) w'^ 

+ ((363238. + 466.224/34 + 49.8562/35 + 7.86612/36 + 1104.3c5,i + 252.84c6,i + 88.4938c7,i) In^ w) w'^ 
+ ((-441763. - 291.503^4 - 34.961/35 - 1215.29c54 - 314.645c64) In^ w) w'^ 
+ ((328765. + 103.587/34 + 699.21c5,i) In'^ w - 130720. In^ w) w'^ 

+ (19838.1 In^ w - 2107.52 In^ w + 99.775 In^ w) w'^ (A.IO) 



-9.125^5 + 0.049/34/35) w-^ 
0.056/39) 

36.489/34 + 0.026/3| - 2.506/35 + 0.033/34/35) w'^ 



11.539 + 14.642/34 - 0.132/31 
0.931/36 - 0.030/37 + 0.198/38 
2124.59 

0.089^6 + 1.336/37 - 0.321/38) w^^ 
(-1406.38 + 37.077/34 - 0.1/3| - 3.784/35 



-680.52 - 60.419^4 + 0.206/3^ 



0.122/36 - 0.803/37 4 
1.535^5 - 2.220 X 10"i%^5) 



0.226/38) In w) 



10.619/36 - 2.373/37 + 2.220 x 10~^'^l3g) w''^ 

(3693.97 + 13.328^4 - 0.117^| - 0.599/35 - 10.012^6 + 2.466^7) In w) w'^ 
(-1841.01 + 11.774/34 - 0.379/35 + 2.497/36 - 0.702/37) In^ w) w"'' 



-15549.1 + 19.182^4 - 3.553 x 10-^^/3^ + II8.7OI/35) 



-1.776 X 10"^'^/34/35 - 26.773/36 - 2.842 x 10-"/37 
(14656. - 22.760^4 - 112.864/35 + 25.512/36) In w) w"^ 
(-3525. + 3.973/^4 + 48.266^5 - 12.021^6) In^ w) w'^ 
(-41.587 + 1.077/34 - 7.103/35 + 1.998^6) In^ w) 
-34138. + 645.226/34 - 

-,-14, 



1.776 x 10"^^/39) 



0.293/3| 



85.394/35 - 0.082/34/35 



4.547 X io-i3« ^--6 



/36j w 

5.684 X IQ-'^'lSj - 3.553 x lO^'^^g - 44.011c5,i + I.975/34C54 - 0.556^5C5,i) 
(53774.4 - 1590.41/34 + 354.169/35) In w + (-23763.2 + 743.3/34 - 168.446^5) In^ w) 
(5321.51 - 195.91/34 + 48.963/35) In^ u; + (-439.615 + 19.731/34 - 5.549/35) In'^ w) u;"^ 



'/36) 
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133168. - 491.03/34 - 0.510^| + 13.826/35 - 0.099/34/35 + 9.095 x 10' 
5.684 X IQ-^^Py + 394.93C54 - 1.317/34C5,i + 124.437c6,i - 0.889/34C6,i) 
(-438674. + 1890.14^4 + 1.171^| - IIO6.IIC54 + 7.901/34C5,i) In w) w"'^ 
(479798. - 2917.82^4) In^ w + (-147748. + 914.619/34) In^ w) w'"^ 
(28316.3 - 187.245/34) In'' w + (-2209.73 + 15.785/34) In'^ u)) w"'^ 

-^^I3l - 11.56^5 - 3.553 x IQ-^^P^P^ - 3.333^6) 

15/ 



-68009.4 - 131.957^4 + 2.842 x 10" 
-2.842 X 10""/37 - 1.421 x IQ-^^Pa - 7.105 x 10" 



^/39 w- 



-231.071C54 - 21.060c6,i - 37.491c7,i) 

(281901. + 257.166/34 + 44.434/35 + 779.652c5,i + 399.906c6,i) In w) w'^ 
(-515266. - 263.312/34 - 1777.36c5,i) In^ w + 454897. In^ w) vj-^ 



-106849. In'' w + 17222.1 In^ w - 1183.57 In^ w) w'^ 

-60068.5 - 439.01/34 - 3.553 x lO'^'^Pl - 33.859/35 - 3.553 x 10" 



-3.402/36 - O.777/37I.42I X IQ-^'^Pa + 1.776 x 10-^^/39 - 546.601c5,i) w 



/34/35) w 
9 



-107.297c6a - 16.309c7,i 
-107.297C64 - 16.309c7,i 



10.485c8,i - 8.882 x 10^1^/3505. 
10.485c8,i)w-^ 



1 w 
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+ ((535556. + 687.396/34 + 73.508/35 + 11.598/36) \nw) w'^ 
+ (1628.16c5,i + 372.785c6,i + 130.475c7,i) Inw) w'^ 

+ ((-976999. - 644.685/34 - 77.318;35 - 2687.73c,5,i - 695.865c6,i) In^ w) w'^ 

+ ((969457. + 305.455/34 + 2O6I.82C54) In^ w - 481830. In"^ w + 87747.5 In^ w - 10875.6 In^ w + 588.427 In^ w) uj-^ 
+ (43141.6 + 479.855/34 + 0.118/3^ + 8.381^5 + 0.099/34/35 + 21.868^6 + 1.844^7) w''^^ 
+ (0.124^8 + 0.056/39 + cioa + 426.754c5,i - 0.658/34C5,i) w'^^ 
+ (0.556/35C5,i + 3.919C6.1 + 0.889/34C6,i + 53.80c7,i + 10.485c84) 

+ ((-449500. - 1271.94^4 - 0.953^| - 509.968/35 - 42.321/36 - 4.255/37 - 1.016^8) In w) ur^^ 

+ ((-1301.71C5.1 - 11.852/34C5,i - 996.586c6,i - 186.115c7,i - 25.284c8,i - 16.C9,i) \nw) w'^^ 

+ ((1.298 X 10^ + 5592.95/34 + 456.071/35 + 50.903/36) In^ w) w"^" 

+ ((8.428^7 + 7894.87C54 + 1613.87c6,i + 334.31c7,i + 113.778c8,i) In^ w) w~^° 

+ ((-2.723 X 10^ - 3004.75/34 - 328.052/35 - 41.953) lii^ w) 

+ ((/36 - 8050.1c5,i - 1907.85c6,i - 471.96707^) In^ w) w'^^ 

+ ((2.352 X 10^ + 1350.17/34 + 139.842/35 + 6104.22c5,i + 1258.58c64) In'' w) w'^^ 
+ ((-1.284 X 10^ - 331.478/34 - 2237.47c5,i) In^ w + 383853. In^ w) 

+ (-49091.1n^w + 4392.42 In^u;- 177.377 ln^w)w"i°. (A.ll) 
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